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Exact solutions are presented for a doubly-periodic array of steadily moving bubbles 
in a Hele-Shaw cell when surface tension is neglected. It is assumed that the bubbles 
either are symmetrical with respect to the channel centreline or have fore-and-aft 
symmetry, or both, so that the relevant flow domain can be reduced to a simply 
connected region. By using conformal mapping techniques, a general solution with 
any number of bubbles per unit cell is obtained in integral form. Several examples 
are given, including solutions for multi-file arrays of bubbles in the channel geometry 
and doubly-periodic solutions in an unbounded cell. 
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1. Introduction 

Since the seminal papers by Saffman & Taylor (1958) and Taylor & Saffman (1959) 
the motion of bubbles in a Hele-Shaw cell, where the fluids are confined between two 
closely spaced parallel plates, has attracted a great deal of attention. Because the 
problem is amenable to analytic treatment, particularly in the case when surface 
tension effects are neglected, many exact solutions have been found for both steady 
and time-dependent flows. A comprehensive list of references, up to 1998, on Hele- 
Shaw flows was compiled by S. D. Howison (people.maths.ox.ac.uk/howison/Hclc- 
Shaw). Hele-Shaw flows have also recently been recognized as having deep con- 
nections with other areas of mathematical physics (Mineev-Weinstein et at 2000, 
Gustafsson & Vasil'ev 2006). Exact solutions for Hele-Shaw flows are thus of inter- 
est not only for their eventual physical significance but also on their own right for 
they are solutions to a difficult free-boundary problem. Some known exact solutions 
for the motion of bubbles and fingers in a Hele-Shaw cell were briefly reviewed by 
Vasconcelos (2001). More recently, other exact solutions have been found, such as 
time-dependent solutions for Hele-Shaw flows around a wedge (Cummings 1999, 
Richardson 2001b, c, Vasconcelos 2007) and solutions for multiple steady bubbles in 
a Hele-Shaw cell with no assumed symmetry obtained by Crowdy (2009a, b). 

Exact solutions for a periodic stream of bubbles in a Hele-Shaw cell have also 
been found. The first of such solutions was obtained by Burgess & Tanveer (1991) 
for a stream of bubbles in the channel geometry with one bubble per unit cell. In 
the limit that the period is taken to infinity their solution recovers the solution 
originally given by Taylor & Saffman (1959) for a single bubble in a channel. Exact 
solutions for a stream of bubbles in an unbounded Hele-Shaw cell were subsequently 
found by Vasconcelos (1993). The Burgess- Tanveer solution was later generalized 
by Vasconcelos (1994) for a wider class of solutions with an arbitrary number of 
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symmetrical bubbles per unit cell. Here again taking the period to infinity yields 
solutions with an arbitrary but finite number of symmetrical bubbles in a Hele- 
Shaw channel (Vasconcelos 200f). It is worth mentioning that a periodic solution 
in a Hele-Shaw channel can be extended to the entire plane by successive reflections 
at the channel walls thus yielding a doubly-periodic solution in an unbounded cell. 
In this context, exact solutions for doubly-periodic flows in a Hele-Shaw cell are of 
interest not only from a mathematical viewpoint, for they correspond to a more 
general situation, but also from a physical perspective because in certain cases they 
can be realized as a stream of bubbles in a Hele-Shaw channel and hence may 
(in principle) be susceptible to experimental investigation. A brief experimental 
study of a stream of bubbles rising in an inclined Hele-Shaw cell was reported by 
Maxworthy (1986). 

In the present paper, we report a rather general class of exact solutions for a 
doubly-periodic array of steadily moving bubbles in a Hele-Shaw cell. Our generic 
solutions are valid in an unbounded Hele-Shaw cell but there is a large subclass of 
solutions that correspond to a stream of bubbles in a Hele-Shaw channel. To render 
the problem analytically tractable, we neglect surface tension effects and suppose 
that the bubbles either are symmetrical with respect to the channel centreline or 
have fore-and-aft symmetry (or both), so that the flow domain can be reduced to 
a simply-connected unit cell. By using conformal mapping techniques, a general 
solution for any number of bubbles per unit cell is obtained in integral form. Our 
solution represents the most general periodic solution for a stream of bubbles in a 
Hele-Shaw cell, in the sense that all previously known periodic solutions are partic- 
ular cases of the solutions reported here. Furthermore, our solutions include several 
novel configurations that were not obtained before, such as the case of a staggered 
two-file array of bubbles in a channel that somewhat resembles the "zipper" -like 
flows of red cells in capillaries (Sugihara-Seki & Fu 2005). Also presented are solu- 
tions with "mixed symmetry" in the channel geometry where within a given bubble 
configuration some bubbles have fore-and-aft symmetry and others are symmetrical 
about the channel centreline. Examples are also given of solutions that cannot be 
restricted to the channel geometry and hence must necessarily be considered in an 
unbounded Hele-Shaw cell. 



2. Mathematical formulation 

We consider the problem of a doubly-periodic array of bubbles moving with a 
constant velocity U along the cc-direction in an unbounded Hele-Shaw cell. For 
convenience, we shall work in a reference frame moving with the bubbles in which 
the flow is steady and the bubbles stationary. We assume that the solutions have 
streamwise and spanwise periods of 2L and 2a, respectively, and so we reduce the 
problem to a rectangular unit cell, as shown in Fig. [T] We place the origin of our 
system of coordinates at the centre of the unit cell, so that its upper and lower edges 
are at y — ±<z and the left and right edges at x = ±L. We shall assume furthermore 
that the bubble configuration is symmetrical with respect to reflections about the 
following four axes: x = 0, x = L, y = 0, and y — a. With this assumption we only 
need to consider the problem in a reduced unit cell corresponding to one quarter, 



say, the upper-right quarter, of the original unit cell; see Fig. 2(a) 
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Figure 1. The unit cell for a doubly-periodic array of bubbles in a Hele-Shaw cell. 



Now let V denote the domain occupied by the fluid within the reduced unit cell. 
For simplicity of notation, we shall denote by C a generic bubble-fluid interface. As 
is usual for Hele-Shaw flows, we introduce the complex potential W(z) = <f> + iip, 
where z = x + iy, 4>{x,y) is the velocity potential, and ip(x,y) is the associated 
stream function. We recall that the velocity potential (in the moving frame) is 
given by (f>(x,y) = — (b 2 pjYl\i) — Ux, where b is the gap between the cell plates, /i 
is the viscosity, and p is the fluid pressure. The complex potential W{z) must be 
analytic in the fluid domain T> and satisfy the appropriate boundary conditions, as 
we now discuss. 

The pressure inside a bubble is assumed to be constant and surface tension 
effects are neglected. The complex potential W must then satisfy the following 
condition on a fluid-bubble interface C 



W = -Ux + C on C, (2.1) 

where C is a complex- valued constant. (The value of the constant C depends on the 
specific bubble one considers and is related to the position of the bubble centroid.) 
The complex potential must also satisfy specific conditions on the boundaries of 
the unit cell, as follows. Given that the upper and lower edges of our unit cell are 
streamlines of the flow and that the left and right edges are equipotentials, we must 
then have 

lmW = on y = 0, (2.2) 

JmW={V-U)a on y = a, (2.3) 

ReVK = on x = 0, (2.4) 

ReW^ = —VL on x = L, (2.5) 

where Im and Re denote real and imaginary parts, respectively. Here V is the 
average velocity across the unit cell in the x-direction in the lab frame and V is 
a positive constant whose physical meaning will become clear shortly. An explicit 
expression for V in terms of L, V and the total area J occupied by the bubbles 
(within a unit cell) can be obtained. We refer the reader to Burgess and Tanveer 
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Figure 2. The reduced unit cell for a doubly-periodic array of bubbles: (a) the 2-plane, 
(b) the VK-plane, (c) the VF-plane, and (d) the ("-plane. 



(1991) for the details and simply quote the results here: 



v = u-v 



UJ 
aL 



(2.6) 



The complex potential W(z) defined above can be seen as a conformal mapping 
from the fluid region in the z-plane onto the corresponding flow domain in the 
W^-plane. In view of Eqs. (|2.1[) - ([2.5[) . one readily sees that the flow domain in the 
W^-plane corresponds to a rectangle, — VL < <j) < 0, — (U — V)a < tp < 0, as 



shown in Fig. 2(b) Bubbles placed along the left and right edges of the unit cell 
in the z-plane are mapped onto horizontal slits in the VF-plane, whereas bubbles 
on the streamlines y — and y — a are mapped onto segments along the lines 
if) = —a(U — V) and ip = 0, respectively; see Fig. 2(b) Note that the horizontal 



slits on Fig. 2(b) correspond to the respective halves of the bubbles on the left 
and right edges that are within the reduced unit cell and hence such slits always 
penetrate into the rectangular domain in the W-plane. 

As discussed in detail elsewhere (Tian & Vasconcelos 1993, Vasconcelos 2001), 
steady Hele-Shaw flows have an interesting rotation invariance in the following 
sense: if a curve C is a solution for a bubble moving with constant velocity U along 
the x-direction in an unbounded Hclc-Shaw cell, then the curve C obtained from a 
rotation of C about the origin by an angle a is also a solution with the same velocity 
U. Alternatively, one can view the rotated solution as one in which the bubble itself 
is kept fixed while its velocity is rotated by an angle a. In view of the rectangular 
geometry of our unit cell, it is particularly relevant for us here to consider a rotation 
by 90°, where the bubbles move with the same velocity U as in the original solution 
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but now in the y-direction. In this case, the region occupied by the fluid in the 
rotated solution is exactly the same as that for the original problem, the main 
difference being that the upper and lower edges of the unit cell which are streamlines 
in the original solution become equipotentials in the rotated problem, whereas the 
left and right edges which were originally equipotentials become streamlines. 

As shown by Vasconcelos (2001), the complex potential W(z) for the rotated 
problem (in the corresponding moving frame) is given by 

W{z)=\[W{z) + Uz]. (2.7) 

From (|2.ip - ([2.5[) and (I2.7[) it then follows that the rotated complex potential W(z) 
satisfies the following boundary conditions: 

W = -Uy + iC on C. (2.8) 

ReW = on y = 0, (2.9) 

ReW = -Va on y = a, (2.10) 

Imf = on i = 0, (2.11) 

lmW=(U-V)L onx = L, (2.12) 

From (|2.12[) one now sees that the parameter V corresponds to the average fluid 
velocity (in the lab frame) in the y-direction of the rotated flow. The flow domain 
for the rotated problem in the VF-plane corresponds to a rectangle, — Va < <f> < 0, 
UL — (f>o < ip < 0, with horizontal slits corresponding to the bubbles placed on 



y = and y = a, as shown in Fig. 2(c) [For the same reasons as discussed before, 
these horizontal slits always penetrate into the bounded domain in VF-planc.] 

Consider now the conformal mapping z(£) from the upper half-£-plane onto 
the fluid domain in the z-plane, so that the boundary dT> of the fluid domain in 
the z-plane is mapped onto the real axis of the £-plane, as shown in Fig. |2(d)| In 
particular, a bubble C will be the image under z(£) of a corresponding interval I 
on the real-£-axis (see below). For dcfinitcncss, we also choose to map the point 
£ = oo to an arbitrary point on the upper edge between the left-upper corner and 
the next bubble; see Fig. 2. In the most general configuration, one has one bubble 
on each corner and an arbitrary number of bubbles on each edge of the cell. To 
keep track of all the bubbles, let us denote by Cf, i — 1,...,4, the bubble on the 
ith corner, with the corners numbered anti-clockwise starting from the upper left 
corner. Under the map z(£) a corner bubble Cf will be image of the interval If on 
the real-£-axis defined by 

J (ai,0i), i = 2,4 , . 

1 ~ I (ft, «,), i = l,3. ( j 

(The reason for distinguishing between even and odd corners is one of convenience, 
for with the above choice one has that in the VK-plane the corners are the images 
of the points ( — cti, whereas in the VF-plane the corners corresponds to £ = ft; 
see Fig. O) Note that if we want a configuration with no bubble in the i-th corner, 
for a given i, we must simply set on — ft. 

Similarly, we shall denote the bubbles on the edges of the unit cell by Cf , where 
the superscript refers to the corresponding side, with a = l,r,u,d, indicating left, 
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right, up and down, respectively, and i — l,...,n Q , where n a is the number of 
bubbles on the respective edge. In the £-plane a given bubble Cf will correspond to 
an interval If on the real axis defined by 

Zr = K 4 _ l7 ^), i = l,...,n a . (2.14) 

Owing to the three degrees of freedom allowed by the Riemann's mapping theorem, 
we could arbitrarily fix the values of three parameters, chosen among a*, ft or vf. 
At this stage, however, it is best to consider all such quantities as free parameters, 
deferring to the examples the moment when we shall need to fix the appropriate 
parameters. 

With some abuse of notation let us now write 

W(0 = W(z(0) and W(0 = W(z(0), 

so that the functions W(C) and W{Q can be viewed as the conformal mappings from 
the upper half-£-plane onto the flow domains in the W- and VF-planes, respectively. 
It then follows from (|2.7p that the mapping z{Q can be written as 

z(Q = ~[w({)+iW(C)]. (2.15) 

Once the mappings W(() and W(C) are known, the bubble interface Cf can then 
be readily computed from (|2.15p by setting £ = s, for s € If. 

Summarizing our procedure so far, we have reduced our original free-boundary 
problem to the much simpler problem of obtaining two conformal mappings, namely, 
W(C) and W((), from the upper half-£-plane onto respective rectangular domains in 
the W- and PF-planes. As we will see shortly, such mappings can be easily obtained 
from the Schwarz-Christofiel formula (Carrier et al. 1983). It should also be noted 
that in previous solutions that use conformal mapping (see, e.g. Tanveer 1987), 
once the mapping W(C) from the chosen domain in the C-plane to the flow domain 
in the M^-plane is known, then the mapping z = /(C) is constructed explicitly so as 
to satisfy the appropriate boundary conditions. Our method has the advantage of 
making this step rather straighforward by identifying the function W(z) = i\W(z) + 
Uz] with the complex potential for the rotated problem, so that solving for the 
mapping W(C,) then completes the solution. Before presenting our general solutions, 
however, it is worthwhile to point out some special features of the solutions for the 
the case U — 2V. 

Solutions with U = 2V are special in the sense that solutions for any U > V can 
be generated by a proper rescaling of the solutions for the former case (Vasconcelos 
2001). To show this, we first note that from inspection of the flow domain in the 
W-plane [Fig. 2(b)| it follows that the mapping Wu(C)i f° r given a and any U > V, 



can be obtained from the corresponding mapping W2v{C), f° r U = 2V, by the 
relation 

Wu(0 = ^^W 2V (0- (2-16) 
In particular, the velocity V [see Fig. [2(b)] for both cases are related by 

U _ v ~ 

Vu = V 2V . (2.17) 
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Similarly, inspection of the flow domain in the M^-plane [see Fig. 2(c) reveals that 

Wu(()=W 2V (0, (2.18) 

with 



(U - Vu)Lu = (2V - V 2V )L 2V . 



(2.19) 



Inserting (|2.17l) into (12 . 19[) yields the half-period Ljj in terms of the corresponding 
value L 2 v for the case U = 2V: 



V{2V - V 2V ) 
UV -{U - V)V 2V 



(2.20) 



It then follows from (|2.15|) . (|2.16p and (|2.18l) that for any U > V the mapping 
Zu(C) can be written in terms of the solutions for U — 2V as 



MO 



1 

u 



1 )W 2 v(0+iW 2V (0 



which can be more conveniently rewritten as 

1 



2V 



MO = -77f7 (1 + P)W 2 v(C) + i(l - P)W 2V (C) 



where 



P = l- 



2V 



(2.21) 



(2.22) 



(2.23) 



with p G (—1,1). This result thus shows that solutions for any U > V can be 
obtained by an appropriate rescaling of the solution with U = 2V (Vasconcelos 
2001). In view of this property, we shall only consider the case U = 2V in the 
remainder of the paper. Without loss of generality we will also set V = 1 and a = 1 
henceforth. 



3. General solutions 



As indicated in Fig. 2(b) the flow domain in the M^-plane is a rectangle with n; 
(n r ) horizontal slits corresponding to the bubbles on the left (right) side of the unit 
cell. This domain can be viewed as a degenerate polygon and so the corresponding 
mapping W(£) can be obtained from the Schwarz-Christoffcl formula (Carrier et 
al. 1983). One finds that W{Q is determined by 



dW 



= K 



4 

n 

;=i 



n(c-7Dn(c-7p 

i=l 3=1 



n(c-oi) n (c-^) n (c-^) 



3 = 1 



k=l 



1/2 ' 



(3.1) 



where K and 7" are real-valued constants to be determined later. In the VF-plane 
the flow domain is of the same form as that in the W-plane, namely, a rectangle 
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with rid (n-u) horizontal slits emanating from the right (left) side of the rectangle. 
By the same reasoning as before, we obtain that the map W(C) is determined by 



dW 
~dZ 



-iK 



n(c--tf) u(c-7f) 

i=l 3 = 1 



4 2n d 2n„ 

n (c - Pi) n (c - v*) n (c - k) 

i=l j=l k=l 



1 1/2 



(3.2) 



where again K and 7? are real-valued constants to be determined. [Note that in 
(|3.f I) and (|3.2j) a given constant 7" is the pre-image of the point on the bubble Cf 
that lies furthest away from the respective edge; see Fig. |2J] 

For calculation purposes it is convenient to expand the numerator in (|3.1D and 
(|3.2I) . and so we write 



dW 
~d( 



dW 
~d( 



K- 



ni+n r 

E «,C 3 



2n; 2n r 

n(c-«i)n (c-^) n (c-^) 
3=1 k=i 



4 

n 

i=i 



1/2 



(3.3) 



n d +n u 



E &K J 



-iK- 



3=0 



4 2n d 2n„ 

n(c-ft)n (c-^) n (c 

i=l j=l fc=l 



1/2 



(3.4) 



1. 



where the Oj's and fej's are real- valued coefficients, with a ni + nr 
[The ay's and fcj's could of course be expressed in terms of the 7j a 's but it is more 
convenient to work directly with the mappings (|3.3p and p. 41) .] For a given set of 
parameters a^, Pi and i/f, the coefficients aj and bj arc determined as follows. 
First, let us define the quantities 7° as 



I a = 

± %3 



Pdt 



' 4 2n, 

n \t-a k \ n 1* 

,k=l m=l 



2n r 



n 



1/2 ' 



(3.5) 



for a — l,r, with i = 1, n a and j = 0, 1, n; + n r . We then note from Fig. [5] 
that the mapping functions W(Q must satisfy the condition W(^2i-i) = WO^fi): 
for a = : Z, r, which in view of (I3.3[) and (I3.5[) yield 



nj+n r 



Z(" aj = 0, for a = l,r and i = 1, 2, n a . 

3=0 



(3.6) 



Now recalling that a„ l+ „ r = 1, we see that p.6[) gives a system of ri; + n r lin- 
ear equations for the n; + n r unknown coefficients ay. Similarly, if we define the 
quantities 



l a — 



Pdt 



' 4 2n d 2n„ 

n n \t-<\ n i*-^ 

1 m— 1 n — 1 



-,1/2 



(3.7) 
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for a — d,u, with i = 1, n a and j — 0, 1, rid + and use the fact that W(C,) 
must satisfy the condition T^(^2i-i) — ^{ v 2i)i f° r a = d,u, we obtain 



for a — u,d and £ = 1, 2, n a , 



(3.8) 



which represent a system of n u + rid linear equations for the n u + rid coefficients bj . 
Thus, once the set of parameters a^, fli and vf are given we can readily obtain the 
constants aj and bj by solving the linear equations (13.6[) and (I3.8p . 

The remaining constants K and K are determined by the boundary conditions 
(|2.3p and (|2.10p . respectively. To see this, first note that K can be calculated from 
the requirement that as we go from point B to point F in Fig. |2(d)| the complex 
potential W in Fig. 2(b) must vary by i (recall that we have made U — 2, V = 1, 
and a = 1), and so we write 



which yields 



K~ 



[Wf B = W(F) - W{B) 



ni+n T 

Y] djP dt 

3=0 



(3.9) 



4 2m 

u(t-»i) n (t 

i=i j=i 



2n r 

v)) n (* 

fe=i 



1/2- 



(3.10) 



The constant if follows from similar requirement: 



IF 



1F(G) -W(A) = 1, 



which yields 

















3=0 





4 2nd 2n„ 

n it-Ai n |*-^| n i*-^ 



i=i 



fc=i 



1/2 



(3.11) 



(3.12) 



We have thus seen that the generic solution described in (|3.3[) and (|3.4j) is 
completely specified by prescribing the set of free parameters a^, and vf, from 
which all other constants appearing in the solution can be calculated. (Recall that on 
account of the three degrees of freedom allowed by the Riemann mapping theorem 
we can fix any three of the above free parameters.) We note furthermore that 
all physical parameters of the solution can be calculated once the mathematical 
parameters are given. For example, the streamwise half-period L can be calculated 
from the requirement that as we go from point F to point M in Fig. 2(a)| the 
mapping z(£) must change by L: 



M 



Am - Ai 



L 



(3.13) 
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Figure 3. Single-file stream of bubbles in a Hele-Shaw channel corresponding to the 
Burgess- Tanveer solution. Here the parameters are 0:2 = —0.6, fa = 0.7, and 

«3 = fa = 0.9. 
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Figure 4. Solutions with (a) centreline symmetry and (b) fore-and-aft symmetry. The 
parameters in (a) are «2 = —0.9, fa = —0.5, uf = —0.49, z/f = 0.4, and Q3 = /?3 = 0.5; 
whereas in (b) we have v\ = —0.9, v\ = —0.6, Q2 = —0.59, fa = —0.2, and Q3 = fa — 0.5. 

which in view of (|2.15[) reads 



[W 



M 



w 



M 



2L, 



(3.14) 



so that from the knowledge of W(() and W(() one readily obtains L. Similarly, 
the bubble centroids and areas can in principle be obtained by carrying out the 
appropriate calculations, but we shall not go into these details here. In the next 
section we will give specific examples of the generic solutions described above. 



4. Examples 

In this section we shall consider some particular cases of the generic solution given 
in the preceding section. We start discussing the situation where we can confine the 
solutions to a Hele-Shaw channel. After this, we will give an example of solutions 
that cannot be realized within the channel geometry and hence must be considered 
in an unbounded Hele-Shaw cell. 
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Figure 5. Solution for a two-file array of bubbles in a channel. Here the parameters are 
u[ = -0.9999, v\ = -0.3, a 2 = #2 = -0.2, a 3 = /3 3 = 0.3, v{ = 0.3001, i/J = 0.7. 
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Figure 6. Solutions for a multi-file array of bubbles with mixed symmetry. The parameters 
are v{ = -0.97, v\ = -0.7, a 2 = /3 2 = -ot 3 = -fa = -0.5, v{ = -0.4, i4 = 0.3, v\ = 0.55, 
and v 2 = 0.9. 



(a) Periodic bubbles in a channel 

Here we consider the particular and experimentally relevant case of a periodic 
array of bubbles moving in a Hele-Shaw channel. We suppose the channel walls are 
at y = ±1, so that y = corresponds to the channel centreline. For this case we set 
in the generic solution given in (|3.3p and (|3.4I) the following conditions 



a\ = Pi = -1, 



Q4 = 04 = 1, 



n u = 0, 



(4.1) 



since there can be no bubble attached to the channel walls. Thus, for the case 
of a periodic array of bubbles in the channel geometry the solution is completely 
specified by prescribing the following set of free parameters: {0^)^=2,3, {/3i}i=2,3, 
and {^ a }i=i,...,2n Q , with a = l,r,d. (Recall that from the third degree of freedom 
allowed by the Riemann mapping theorem one of these parameters can be kept 
fixed.) We also note that the generic solutions for multiple bubbles in a Hele-Shaw 
channel presented by Vasconcelos (2001) can be obtained as special case of our 
periodic solution above in the limit L — > 00. This can be accomplished by simply 
setting Pi = P2 = — 1 and P3 — P4 = 1 in Fig. 2(d) In what follows we shall discuss 



some specific examples of periodic solutions in a Hele-Shaw channel. 

We start with the case of a single bubble per unit cell. In the notation of Fig. [1] 
this is achieved by placing a bubble at the lower left corner of the reduced unit cell 
and no other bubble anywhere. In this case we reproduce the exact solution for an 
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Figure 7. Solutions with mixed symmetry including a symmetrical bubble at the centre 
of the unit cell. The parameters are v\ — —0.8, v\ — —0.5, a-z = —0.48, fa = —0.1, 
v{ = -0.08, v$ = 0.3, i/f = 0.32, vf = 0.6, and a 3 = fa = 0.8. 



infinite stream of bubbles in a Hele-Shaw cell originally obtained by Burgess and 
Tanveer (1991) using a rather different method. For completeness we show in Fig. [3] 
an example of the Burgess- Tanveer solution. Note, in particular, that in this case 
the bubble is placed at the centre of the unit cell and therefore has both centreline 
and fore-and-aft symmetry. 

As the next example, we consider the situation in which there are bubbles only 
either at the lower edge or at the left edge of our reduced unit cell. Examples of 
these two cases are shown in Figs. 4(a) and |4(b)| In such cases we recover the class of 
periodic solutions with centreline symmetry or fore-and-aft symmetry, respectively, 
obtained earlier by Vasconcelos (1994). It should be noted, however, that in the 
work by Vasconcelos (1994) no specific example was given, only the general method 
of solution was presented. In this sense, this is the first time that explicit solutions 
for such symmetrical configurations have been calculated. 

Next we note that novel solutions can be generated if we place bubbles simulta- 
neously at the left and at the right edges of the reduced unit cell (and nowhere else). 
For example, in Fig. [5] we show a solution with two bubbles per unit cell, where 
each bubble has fore-and-aft symmetry but they lie in different equipotentials of the 
flow. Seen in the laboratory frame, this solution corresponds to a staggered two-file 
array of bubbles moving steadily down the channel. It is interesting to notice that 
similar zipper-like arrangements are observed in the flow of red cells in capillaries 
(Sugihara-Seki k Fu 2005). 

Another type of new solutions includes the case of "mixed symmetry" where 
some bubbles have centreline symmetry and some others have fore-and-aft sym- 
metry. An example of such case is the multi-file array of bubbles shown in Fig. [6j 
where there is one file of bubbles symmetrical about the channel centreline, with the 
other files of bubbles having fore-and-aft symmetry. Another solution with mixed 
symmetry is shown in Fig. [7] where in this case there is a symmetrical bubble at 
the centre of the unit cell. 
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Figure 8. Solution for a doubly-periodic array of bubbles in an unbounded Hele-Shaw 
cell. Here the parameters are /3i = —1.3, ai = — 1, v\ — —0.9, v\ = —0.6, a 2 = —0.5, 
/3 2 = -0.4, uf = -0.3, ui = 0.3, /3 3 = 0.32, a 3 = 0.5, v\ = 0.55, u r 2 = 0.8, a 4 = 1, 
/3 4 = 1.1, i/f = 1.12, and u% = 1.5. 

(b) Solutions in an unbounded cell 

Here we briefly discuss the case in which the periodic array of bubbles cannot 
be physically realized within the channel geometry and hence must be considered 
in an unbounded Hele-Shaw cell. We begin by noting that if there are bubbles 
placed on the upper and lower edges of our reduced unit cell, then the solutions 
cannot be reduced to the channel geometry as before, and so one must considered 
an unbounded Hele-Shaw cell. An example of such case is given in Fig. \8\ It is 
worth mentioning that this type of solution represents the most general solution for 
a periodic array of (symmetrical) bubbles in a Hele-Shaw cell, in the sense that the 
other types of solution can be obtained by either 'decimating' or 'adding' bubbles 
along any of the edges of the reduced unit cell. 



5. Conclusions and Discussion 

We have presented an exact solution for a doubly-periodic array of steadily moving 
bubbles in a Hele-Shaw cell with an arbitrary number of bubbles per unit cell. 
Our solution represents the most general periodic solution for bubbles in a Hele- 
Shaw cell known to date, in the sense that the previously known solutions for 
a stream of bubbles in a Hele-Shaw channel are particular cases of the solutions 
reported here. In addition, our solution includes novel cases of periodic arrays of 
bubbles in a Hele-Shaw channel, such as multi-file flow of symmetrical bubbles and 
configurations with mixed symmetry, where some bubbles have centreline symmetry 
while others have fore-and-aft symmetry. Examples of a doubly-periodic array of 
bubbles in an unbounded cell that cannot be restricted to the channel geometry 
were also discussed. 

In constructing our solutions we assumed that the bubbles either are symmet- 
rical either with respect to the cell centreline or have fore-and-aft symmetry (or 
both) , so that the relevant flow domain can be reduced to a simply-connected unit 
cell. In this context, the existence of periodic solutions in a Hele-Shaw cell with- 
out any symmetry requirement, where the unit cell is no longer simply-connected, 
remains an open and interesting question. Recently, a new methodology to treat 
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Hele-Shaw flows in multiply-connected domains has been developed by Richard- 
son (2001a) and Crowdy (2009a, b). In particular, solutions for a finite number of 
asymmetric steadily-moving bubbles have been found both in an infinite Hele-Shaw 
cell (Crowdy 2009a) and in the channel geometry (Crowdy 2009b). It is therefore 
likely that our solutions can be generalized to include the case of a periodic array of 
asymmetric bubbles. This interesting (and more difficult) problem clearly deserves 
to be studied further. 

Another interesting point to be investigated is the effect of surface tension on 
the shape and velocity of our periodic array of bubbles. From previous works on 
the related problem for a single bubble (Tanveer 1986, 1987; Combescot & Dom- 
bre 1988) as well for the Burgess- Tanveer periodic solution (Burgess & Tanveer 
1991), one expects that the 'degeneracy' displayed by our solutions, meaning that 
fixing the geometrical parameters of the solution does not fix the bubble veloc- 
ity U, will most certainly be removed by surface tension. However, it is not clear 
which particular bubble arrangements, among our large class of solutions, will have 
counterparts in the case of nonzero surface tension. For instance, in the case of a 
single bubble in a Hele-Shaw channel, Tanveer (1987) was unable to find solutions 
without symmetry with respect to channel centreline, even though such nonsym- 
metric solutions exist when surface tension is neglected. In light of these results, it 
appears that periodic solutions with centreline symmetry are more likely to 'sur- 
vive' in the nonzero surface tension scenario. Unfortunately, however, the search 
for solutions with surface tension has to be performed numerically — a task that is 
likely to become increasingly more demanding as the number of bubble grows. 

This work was supported in part by the Brazilian agencies FINEP, CNPq, and FACEPE 
and by the special programs PRONEX and CT-PETRO. 

References 

Burgess, D. & Tanveer, S. 1991 Infinite stream of Hele-Shaw bubbles. Phys. Fluids A 3, 
367-379. 

Carrier, G. F., Krook, M. & Pearson, C. 1983 Functions of a complex variable. Ithaca: 
Hod Books. 

Combescot, R. & Dombre, T. 1988 Selection in the Saffman- Taylor bubble and asymmetric 

finger problem. Phys. Rev. A 38, 2573-2581. 
Crowdy, D. 2009a Multiple steady bubbles in a Hele-Shaw cell. Proc. R. Soc. A 465, 

421-435. 

Crowdy, D. 2009b An assembly of steadily translating bubbles in a Hele-Shaw channel. 

Nonlineanty 22, 51-65. 
Cummings, L. J. 1999 Flow around a wedge of arbitrary angle in a Hele-Shaw cell. Eur. 

J. Appl. Math. 10, 547-560. 
Gustafsson, B. & Vasil'ev, A. 2006 Conformal and potential analysts in Hele-Shaw cells. 

Basel: Birkhauser Verlag. 
Maxworthy, T. 1986 Bubble formation, motion and interaction in a Hele-Shaw cell. J. 

Fluid Mech. 173, 95-114. 
Mineev-Weinstein, M., Wiegmann, P. B. & Zabrodin, A. 2000 Integrable structure of 

interface dynamics. Phys. Rev. Lett. 84, 5106-5109. 
Richardson, S. 2001a Hele-Shaw flows with time-dependent free boundaries involving a 

multiply-connected fluid region. Eur. J. Appl. Math. 12, 571-599. 



Article submitted to Royal Society 



Doubly-periodic array of bubbles in a Hele-Shaw cell 



15 



Richardson, S. 2001b Hele-Shaw flows with free boundaries in a corner or around a wedge 

Part I: Liquid at the vertex. Eur. J. Appl. Math. 12, 665-676. 
Richardson, S. 2001c Hele-Shaw flows with free boundaries in a corner or around a wedge 

Part II: Air at the vertex. Eur. J. Appl. Math. 12, 677-688. 
Saffman, P. G. & Taylor G. I. 1958 The penetration of a fluid into a porous medium or 

Hele-Shaw cell containing a more viscous liquid. Proc. R. Soc. Lond. A 245, 312-320. 
Sugihara-Seki, M. & Fu, B. M. 2005 Blood flow and permeability in microvessels. Fluid 

Dyn. Res. 37, 82-132. 

Tanveer, S. 1986 The effect of surface tension on the shape of a Hele-Shaw cell bubble. 

Phys. Fluids 29, 3537-3548. 
Tanveer, S. 1987 New solutions for steady bubbles in a Hele-Shaw cell. Phys. Fluids 30, 

651-658. 

Taylor, G. I. & Saffman, P.G. 1959 A note on the motion of bubbles in a Hele-Shaw cell 

and porous medium. Q. J. Mech. Appl. Math. 12, 265-279. 
Tian, F. R. & Vasconcelos, G. L. 1994 Rotation invariance for steady Hele-Shaw flows. 

Phys. Fluids A 5, 1863-1865. 
Vasconcelos, G. L. 1993 Exact-solutions for a stream of bubbles in a Hele-Shaw cell. Proc. 

R. Soc. Lond. A 442, 463-468. 
Vasconcelos, G. L. 1994 Multiple bubbles in a Hele-Shaw cell. Phys. Rev. E 50, R3306- 

R3309. 

Vasconcelos, G. L. 2001 Exact solutions for steady bubbles in a Hele-Shaw cell with rect- 
angular geometry. J. Fluid Mech. 444, 175-198. 

Vasconcelos, G. L. 2007 Comment on "Free surface Hele-Shaw flows around an obstacle: 
A random walk simulation". Phys. Rev. E 76, 038301. 



Article submitted to Royal Society 



